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1.  GENERAL 

Adaptive  feedback  and  feedforward  control  techniques  can 
significantly  improve  accuracy  and  dynamic  performance  of  the 
existing  pointing-acquisition-tracking  (PAT)  mechanisms  used  in 
laser  communication  systems  [3,4,5].  Although  such  systems 
utilize  the  most  advanced  materials  and  technologies,  their 
performance  is  adversely  affected  by  the  exposure  to  jitter  and 
natural  transients  (bending  modes)  of  moving  mechanical  parts. 

The  main  motivation  for  this  research  is  in  the  fact  that 
improving  the  accuracy  of  a  laser  beam  steering  system  has  the 
potential  for  replacing  an  expensive,  high  performance  optical 
mirror  by  a  lower  cost  mirror  with  enhanced  performance  of 
control  system.  In  addition,  improvement  of  beam  positioning 
accuracy  can  result  in  the  reduction  of  the  required  power  of 
lasers  used  as  communication  medium  [5]. 

The  study  was  aimed  at  the  application  of  advanced  control 
techniques  for  further  improvement  of  dynamic  performance  of  a 
500Hz  fast  steering  mirror  of  a  pointing-acquisition-tracking 
(PAT)  system  for  lasercom.  First,  a  series  of  laboratory  tests 
was  conducted  leading  to  the  analysis  of  dynamic  characteristics 
of  the  mirror.  Recorded  responses  of  the  azimuth  and  elevation 
positions  of  the  mirror  were  utilized  for  the  development  of  a 
detailed  mathematical  model  of  the  mirror,  featuring  "main 
dynamics"  of  the  elevation  and  azimuth  channels,  bending  modes, 
and  cross  coupling  between  the  channels.  The  model  was 
implemented  in  software  and  its  parameters  were  estimated  by 
matching  simulated  responses  to  experimental  data. 

The  mathematical  and  simulation  models  obtained  herein,  were 
used  for  comprehensive  investigation  of  the  mirror  performance 
and  analysis  of  its  compliance  with  design  specifications.  The 
mathematical  description  provided  the  basis  for  the  modification 
of  control  procedures  resulting  in  the  extension  of  performance 
limits  of  the  mirror  hardware.  The  model-based  simulation  scheme 
was  used  to  define  and  evaluate  the  reference  model,  feedback  and 
feedforward  controllers,  assuring  that  improvements  in  mirror 
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performance  are  weighted  against  magnitudes  of  control  efforts. 

The  model  reference  (MR)  approach  of  adaptive  control  was 
used  for  the  modification  of  the  existing  control  system  of  the 
500Hz  mirror.  The  reference  model,  feedback  and  feedforward 
controllers  were  defined  in  terms  of  transfer  functions  and 
software  procedures.  Application  of  the  MR  approach  allowed  for 
the  cancellation  of  undesirable  components  of  mirror  dynamics 
(bending  modes) .  The  MR  approach  provided  the  basis  for  the 
elimination  of  cross-coupling  between  azimuth  and  elevation 
channels.  Frequency  response  of  the  particular  channels  of  the 
mirror  was  flattened  and  significantly  extended.  The  study 
provided  the  evidence  that  scope  of  applications  of  the  500  Hz 
mirror  could  be  widened  by  improving  its  control  system. 

Effects  of  satellite  jitter  on  the  beam  positioning  system 
are  quite  similar  to  the  effects  of  bending  modes  of  the  optical 
mirror.  A  simulation  study  of  the  application  of  the  MR  approach 
for  the  reduction  of  jitter  effects  on  a  high  precision 
positioning  system  was  conducted. 

Preservation  of  performance  characteristics  of  a  PAT  system 
subjected  to  long  term  effects  of  adverse  operational  conditions 
and  natural  wear  presents  an  important  problem.  A  simulation 
study  of  the  application  of  the  MR  approach  for  improving 
robustness  of  a  high  precision  positioning  system  was  conducted. 

The  experimental  part  of  this  study  was  conducted  on  the  PAT 
testbed  of  Rome  Laboratory.  The  computer  simulation  work 
utilized  software  package  CC  providing  an  analytical  solution  to 
the  system  analysis  and  control  problems  [9]. 
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2.  EXPERIMENTAL  ANALYSIS  OF  A  500Hz  OPTICAL  MIRROR 

Mathematical  model  of  a  500  Hz  beam  steering  mirror, 
including  its  control  circuitry,  is  based  on  the  series  of 
transient  responses  recorded  during  laboratory  testing.  The 
testing  procedure  implied  that  step  voltage  signals  were  applied 
to  the  azimuth  (X)  and  elevation  (Y)  inputs  of  the  mirror  control 
system.  The  responses  of  the  azimuth  (X)  and  elevation  (Y) 
positions  of  the  mirror  were  represented  by  the  output  voltage 
signals  of  built-in  sensors  of  the  mirror  system.  Step  responses 
of  the  mirror  channels  were  recorded  repeatedly  at  different 
magnitudes  of  the  input  signals  in  order  to  assess  the  linearity 
of  mirror  characteristics  and  repeatability  of  results. 

Generation  of  the  test  signals  and  recording  of  testing  results 
were  performed  using  a  PC  equipped  with  software  and  interface  of 
LabView  package.  The  scheme  shown  in  Figure  2-1  illustrates  the 
testing  setup. 

Testing  results  are  given  as  follows: 

Figure  2-2a  shows  the  family  of  responses  of  the  azimuth 
position  to  the  .5,  1.,  1.5,  2.,  3.,  and  3.2  V  signals  applied  to 
the  azimuth  input  (X-to-X  response) .  Figure  2-2b  shows  one  of 
the  step  responses  recorded  at  a  much  higher  sampling  frequency. 
Additional  records  of  transient  responses  of  the  azimuth  channel 
are  given  in  figures  A-1,  A-2 ,  and  A-3  of  the  Appendix. 

Similarly,  Figures  2-3a,b  shows  the  family  of  responses  of 
the  elevation  position  to  the  1.,  1.5,  2.,  2.5,  3.,  and  3.5  V 
x’eference  signals  (Y— to— Y  response)  .  Additional  records  of 
transient  responses  of  the  elevation  channel  are  given  in  figures 
A-4,  A-5,  and  A-6  of  the  Appendix, 

Figures  2-3a,b,c,d,f  represent  crosscoupling  between  the 
azimuth  and  elevation  channels,  providing  an  evidence  that 
crosscoupling  can  be  viewed  as  the  "exchange  of  bending  modes" 
between  channels  with  some  residual  effects.  Additional  records 
of  transient  responses  of  the  elevation  channel  are  given  in 
Figures  A-7  of  the  Appendix. 
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Figure  2-2.  Recorded  step  response  of  the  azimuth  channel 
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Figure  2-4.  Recorded  cross -coupling  interaction  betweer 
the  elevation  and  azimuth  channels 
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3.  MATHEMATICAL  DESCRIPTION  AND  COMPUTER  SIMULATION  OF  THE 

500Hz  MIRROR 

Recorded  step  responses  of  the  azimuth  and  elevation 
channels  of  the  mirror  were  sufficient  for  the  evaluation  of  the 
transfer  functions  of  the  mirror  channels.  This  task  was 
performed  by 

-  interpretation  of  the  recorded  responses  as  a  combination 
of  step  responses  of  "typical”  dynamic  systems, 

-  parameter  estimation  of  "typical"  systems  accountable  for 
the  frequency,  overshoot,  and  settling  time  of  each  component  of 
the  recorded  response,  and 

-  estimation  of  "weights"  of  each  component  accountable  for 
the  overall  shape  and  steady-state  value  of  the  recorded 
response . 

The  mathematical  description  task  utilized  analysis  of  step 
response  of  many  "candidate  models"  provided  by  CC  software,  and 
can  be  viewed  as  an  "intelligent  trail-and-error"  process. 

Analysis  of  the  obtained  responses  indicates  that  X-to-X 
dynamics  exhibits  a  combination  of 

-  dynamics  of  a  first  order  system  with  settling  time  of  .02 
seconds , 

-  dynamics  of  a  second  order  system  with  settling  time  of 
.0005  seconds,  natural  frequency  of  3800  rad/sec,  damping  ratio 
of  .5, 

-  dynamics  of  a  second  order  system  with  settling  time  of 
.0095  seconds,  natural  frequency  of  3500  rad/ sec,  damping  ratio 
of  .03,  and  constant  gain  of  0,  and 

-  dynamics  of  an  integrator,  resulting  in  a  ramp-type  drift 
with  the  slope  of  .3  units/sec. 

It  is  believed  that  the  first  two  components  are  responsible 
for  the  "main"  dynamics  of  the  azimuth  channel,  and  can  be 
defined  as 

G^^i(s)=32/(s+200)  , 
and 

G^jj2(s)=l-444e7/  (s2+3800s+l.  444e7)  , 
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The  third  term  represents  a  long-lasting  high  frequency 
vibration  which  propagates  through  the  mirror  structure  (bending 
mode)  and  does  not  results  in  any  residual  effects.  It  can  be 
defined  as 

Gx^3(s)=183.75s/  (s2+210s+1.225e7)  , 

The  drift  term,  clearly  seen  on  the  recorded  step  response 
curves,  can  be  defined  as 

^XX4  (  ~  ‘  ^  ^  ^ 

Finally,  mathematical  model  of  the  azimuth  control  channel 
of  the  500Hz  mirror  (closed-loop)  is 

^XX  ~*^XX1  (  ®  )  '*'‘^XX2  "’'*^XX3  '''^XX4  ‘ 

The  simulated  step  responses,  given  in  Figures  3-1  and  3-3, 
show  very  good  resemblance  of  the  experimental  results. 

Figure  3-2  features  the  frequency  response  of  the  X-to-X 
channel  which  exhibits  noticeable  nonlinearity  within  the 
operational  frequency  range,  and  a  resonant  peak  at  the  frequency 
of  bending  modes. 

Mathematical  description  and  simulation  of  the  closed-loop 
transfer  function  of  the  elevation  channel,  Gyy(s),  is  very 
consistent  with  Gj(^(s)  and  was  also  quite  successful: 

Gyy  (  S  )  ~^yY1  (  ®  )  ^^YY2  (  ®  )  '^^YY4  (  ®  ^  ' 

where 

Gyyt  (s)=32/  (s+200)  , 

GYY2(s)=l-444e7/  (s2+2800s+l .  444e7 )  , 

Gyy3(s)=183.75s/  (s2+210s+1.225e7)  , 

Gyy4  (^)~* 

Mathematical  description  of  the  crosscoupling  channels  was 
obtained  in  the  form, 

G^y(s)=Gyjj(s)=150.s/  (s2+210s+1.225e7) 
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Figure  3-1.  Step  response  of  the  math. model  of  the  X-X  channel 
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Figure  3- 


3.  Step  response  of  the  math. model  of  the  X-Y  channel 
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4.  APPLICATION  OF  THE  MODEL  REFERENCE  APPROACH  FOR 
PERFORMANCE  MODIFICATION  OF  THE  500Hz  MIRROR 

Availability  of  mathematical  models,  describing  the  "main" 
dynamics  of  the  X-to-X,  Y-to-Y,  X-to-Y,  and  Y-to-X  mirror  control 
channels,  facilitates  the  implementation  of  the  model  reference 
(MR)  control  scheme,  eliminating  undesirable  components  of  mirror 
dynamics,  such  as  bending  modes  and  crosscoupling. 

The  block  diagram  below  illustrates  the  principle  of 
operation  of  the  MR  scheme. 


Blocks  X-to-X  and  Y-to-X  represent  dynamics  of  the  "direct" 
mirror  channel,  and  crosscoupling  between  channels.  Block  "RM" 
(reference  model)  represents  the  "main"  dynamics  of  the 
appropriate  channel,  free  of  bending  modes  and  effects  of 
crosscoupling.  It  could  be  seen  that  error  signal  z(t) 
represents  all  undesirable  motion  components  (including  bending 
modes  and  effects  of  crosscoupling) ,  which  are  extracted  by  the 
error  detector.  Block  "compensator"  is  responsible  for  the 
transformation  of  the  detected  error  signal  into  compensation 
signal  u(t),  injected  in  the  reference  path  of  the  mirror  control 
circuitry. 

Design  of  the  MR  system  is  aimed  at  the  selection  of  such 
compensation  law,  H(s) ,  which  would  "force"  the  output  of  the 
mirror  follow  the  output  of  the  reference  model.  Since  the 
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output  of  the  reference  model  includes  only  the  "main  dynamics" 
terms,  the  compensator  is  expected  to  reduce  all  undesirable 
components  of  the  system  response  such  as  bending  modes  and 
effects  of  crosscoupling.  However,  the  benefits  of  a  MR  system 
do  not  come  without  a  price:  compensation  signal  u(t)  results  in 
the  increased  magnitude  of  the  control  efforts  in  the  mirror 
positioning  system.  The  increased  control  effort  may  or  may  not 
be  carried  out  by  the  existing  mirror  drive  devices,  and  may 
results  in  some  additional  explicit  or  implicit  costs. 

Therefore,  the  feasibility  of  control  effort  is  one  of  the 
important  considerations  during  the  design  of  a  MR  system.  While 
the  magnitude  of  control  effort  is  of  primary  concern,  computer 
simulations  present  the  only  opportunity  for  dealing  with  this 
issue . 

It  can  be  shown  that  the  following  is  the  equivalent 
representation  of  the  above  system,  more  suitable  for  computer 
simulation  analysis  and  design  of  the  mirror  control  circuitry: 


The  model  reference  control  scheme,  considered  above,  was 
developed  for  the  500  Hz  optical  mirror  and  implemented  in  the 
form  of  a  simulation  model.  The  following  first  order  system  was 
selected  as  the  reference  model; 

^hooel(s)=4680/(s+4000) 

It  was  found  that  the  most  suitable  controller  H(s)  can  be 
defined  in  the  form  of  a  phase-lead  circuit, 

H (s)=(l+. 001s) / (1+. 0001s) 
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step  response  of  the  X-X  channel,  after  adding  the 
compensation  scheme,  is  shown  in  Figure  4-1.  It  can  be  seen  that 
the  response  is  virtually  free  from  the  high  frequency  "bending 
mode"  terms.  Frequency  response  of  the  X-X  channel  is  shown  in 
Figure  4-2 .  It  can  be  seen  that  application  of  the  MR  approach 
has  resulted  in  the  enhanced  operational  characteristics  of  the 
500Hz  mirror:  its  bandwidth  is  almost  doubled,  and  the  magnitude 
of  the  resonant  peak,  corresponding  to  bending  modes,  is 
drastically  reduced. 

Recall  that  in  the  original  simulation  tests  the  magnitude 
of  the  input  signal,  applied  to  the  mirror  circuitry,  was  equal 
to  1.0.  Figure  4-3  presents  the  control  effort,  generated  in  the 
model  reference  compensation  scheme:  its  magnitude  is  equal  to 
2.,  i.e.  the  "cost"  of  the  performance  enhancement  is  doubling 
the  magnitude  of  the  control  effort. 

Application  of  the  reference  model 

Ghcx)El(s)=‘7°20/(s+6000) 

results  in  the  appropriate  reduction  of  the  settling  time, 
elimination  of  high  frequency  transient  terms,  and  extends  the 
system  bandwidth  to  1.5KHz  with  the  increase  of  the  control 
effort  by  200%.  (See  Figures  4-4  and  4-5). 

In  theory,  a  model  reference  compensation  scheme  is  capable 
of  the  elimination/reduction  of  all  "undesirable"  components  of 
the  transient  response,  including  the  ones  caused  by  cross- 
coupling. 

A  simulation  study  of  the  decoupling  effect  of  the  model 
reference  scheme  was  conducted.  The  block  diagram  of  the 
equivalent  simulation  setup  is  shown  below: 

Y  ref. 


>  X 
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Figure  4-1 
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Figure  4-2 
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Figure  4-3 
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Figure  4-4. 
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Figure  4-6 


REDUCTIOM  OF  THE  CROSS-COUPLING  EFFECT  IN  THE  NR  SCHENE 


It  could  be  seen  that  reduction  of  cross-coupling  effects 
depends  on  the  frequency  response  of  the  transfer  function  with 
unity  forward  path  and  the  feedback 

Gpg(s)=Gx^(s)H(s)  = 

[32/(s+200)+1.444e7/(s2+3800s+1.444e7)+ 

+83.75S/ (s2+210s+1.225e7)+.3/s] [ (1+. 001s) / (1+. 0001s) ] 

The  appropriate  frequency  response  indicates  that  only  modest 
reduction  of  cross-coupling  effect  (5  db)  can  be  expected. 
However,  this  reduction  can  be  achieved  by  increasing  the  gain  of 
the  compensator  H(s) .  For  example,  Figure  5-6  indicates  that  at 

H(s) =10(1+. 001s) / (1+. 0001s) 

cross-coupling  effect  can  be  reduced  by  20  db  within  the 
frequency  range  of  interest.  According  to  simulation  analysis, 
the  compensator  gain  can  be  increased  up  to  100  without  any 
stability  loss  or  noticeable  increase  of  the  magnitude  of  control 
efforts . 

The  implementation  of  the  MR  control  technology  requires 
that  the  model  reference  and  the  compensator  be  properly  defined 
in  software  or  hardware.  The  following  discrete-time  expressions 
can  be  utilized  for  the  software  implementation. 

Reference  model  via  the  central  difference  scheme: 

Y(s) /X(s) =k/ (s+a) ,  or 
s Y ( s ) +a Y ( s ) =kX ( s ) ,  or 
y' (t) +ay (t)=kx(t) ,  or 

y[  (i+1)  A^]-y[  (i-1)  A^]+2aA*y[iA*]=2kA^x[iA^]  ,  or 
Yi+i-i'i-i+2aA^y.=2kA^Xj,  or 
y,.^i=y,..i-2aA^y,+2kA^x,. , 

where 

y(t),  x(t),  Y(s),  and  X(s)  are  output  and  input  signals  of 
the  reference  model  and  their  Laplace  transforms, 

k  and  a  are  parameters  of  the  reference  model, 
t  -  continuous  time, 

A*  is  time  discretization  step,  and 
y[iA^]=y,=y(t) /t=iA^  and  x[iA*]=x,.=x(t) /t=iAS 
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The  compensator  via  the  central  difference  scheme  (using  the 
same  notations) : 

Y(s) /X(s)=k(l+as) / (1+bs) ,  or 
sbY ( s ) +Y ( s ) =ksaX ( s ) +kX ( s ) ,  or 
by' (t) +y (t) =kax' (t) +kx(t) ,  or 
byi^i-by,..i+2A^y.=kax.^^-kax,..,+2kA^ ,  or 
by,--by..2+2A^y..^=kax.-kaXj.2+2kA^Xj.^ ,  or 

yi=[byi.2-2AV{.i+kax.-kax,..2+2kA^x..^]/b 

Implementation  of  the  MR  control  using  analog  circuitry: 


Y(s)  _  1 

X(s)  RCs+l 


Y(s)  s+l/RC 

X{s)  s+(R+r)/RrC 


CONCLUSION:  Characteristics  of  high  performance  optical 
mirrors  can  be  significantly  improved  by  upgrading  their  existing 
control  systems  using  model  reference  control  schemes: 

-  freguency  response  can  be  extended  and  linearized  within 
the  operational  frequency  range, 

-  high  freguency  components  of  the  transient  response 
(bending  modes)  can  be  eliminated, 

-  the  cross-coupling  effect  can  be  controlled. 

The  potential  of  the  MR  approach  is  limited  only  by  the  allowable 
magnitude  of  control  effort. 


-  first  order  system 


phase-lead  system 
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5.  MATHEMATICAL  DESCRIPTION  AND  COMPUTER  SIMULATION  OF 
A  HIGH  PRECISION  POSITION  CONTROL  SYSTEM 


This  section  of  the  Report  presents  a  simulation  study  on 
various  aspects  of  the  application  of  the  Model  Reference  (MR) 
approach  in  high  precision  positioning  systems.  While  in  the 
previously  conducted  studies,  the  original  positioning  system  was 
defined  by  a  single  dynamic  block  representing  ’’closed-loop” 
system  behavior,  it  is  assumed  herein  that  the  entire  block 
diagram  of  the  original  system  and  its  particular  modules  are 
known.  This  assumption  allows  for  the  explicit  and  accurate 
investigation  of  control  efforts  in  the  original  and  MR  systems, 
system  sensitivity  to  jitter,  effects  of  cross-coupling,  and 
system  sensitivity  to  internal  parameters  of  the  controlled 
plant.  Note,  that  although  the  configuration  of  the  system 
described  below  is  typical  of  a  fine  steering  mirror  control 
system,  the  study  is  concerned  with  general  principles  and  does 
not  target  any  specific  beam  steering  technology. 

Consider  a  high  precision  position  control  system 
represented  by  the  block  diagram  below. 

V 


The  control  loop  of  this  system  is  comprised  of  electronic 
and  electro-mechanical  elements.  The  system  dynamics  includes 
high-frequency  components  of  the  transient  process,  so-called 
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bending  modes.  Presently,  bending  mode  components  are  not 
"covered"  by  the  existing  position  sensor.  They  are  not  present 
in  the  feedback  signal,  and  are  not  compensated  for  by  the 
existing  controller. 

Note  that  double  lines  symbolize  vector  signals,  and  blocks 
represent  transfer  matrices, 

vector  X=[X^,X2]^  represents  the  azimuth  and  elevation 
components  of  the  position  is  the  transpose  symbol)  , 

vector  R=[R^,R2]^  represents  the  azimuth  and  elevation 
components  of  the  position  reference, 

signal  V  is  a  scalor  representing  the  vibration  (jitter) 
affecting  the  system, 

transfer  matrix  G(s)  represents  the  dynamics  of  the 
controlled  process  comprised  of  "direct  channel"  transfer 
functions  giT(s)  and  g22(s)  ,  transfer  functions  g^2(s)=g2i  (s) 
representing  crosscoupling  effects,  and  transfer  functions 
g^3 (s)  =g23 (s)  representing  effects  of  vibration. 


G(s)  = 


gii  (s) 
g2i  (s) 


gi2(s) 

g22(s) 


gi3(s) 

g23(s) 


Transfer  matrix  H(s)  represents  the  combined  dynamics  of  the 
system's  drive  generating  the  azimuth  and  elevation  control 
efforts,  U, ,  and  U2,  and  the  prefilter  placed  in  the  system's 
forward  path: 


H(s) 


hii(s)  0 

0  h22(s) 


Transfer  matrix  B(s)  represents  the  effect  of  "bending 
modes"  on  the  position  output. 


B(s)  = 


b„(s) 

0 


0 

b22(s) 
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Transfer  matrix  F(s)  represents  the  feedback  controllers, 

fii(s)  0 
F(s)=  0  f22(s) 

Transfer  matrix  P(s)  represents  the  prefilters  in  the 
reference  channel, 

Pii(s)  0 
P(s)=  0  P22(s) 

For  the  purpose  of  simulation,  the  following  definitions  of 
the  particular  transfer  functions  were  selected: 
gii  (s)=g22(s)  =  (s+2)  /[s(s+20)  ] 

(s)=100s/  (s2+36s+6002) 
b22(s)=126s/  (s2+112s+800^) 
g^2(s)=2.3s/  (s2+30s+1502) 

g,3(s)=310/  (s2+12 . 3  6S+3  600) +66*s/  ( s^+9  .  Ss+SOO^) 
h^,  (s)  =h22  (s)  =  (s+2  0)  /  (s+40) 
f^^  (s)=f22(s)=1000/  (S+2  0) 

Pii  (s)=P22(s)  =1000/  (s+2  0) 

On  the  basis  of  these  definitions  and  system  configuration, 
the  following  transfer  functions  of  interest  can  be  established: 

-  transfer  function  of  the  existing  system  (reference-to- 
position) : 

GE(s)=p^,(s)Go(s)+b^^(s) 

where 

(s)g„  (s) 

Gq  ( s )  = - 

l+hii(s)g^^(s)f„(s) 

-  transfer  function  of  the  existing  system  (ref erence-to- 
control  effort) : 

p„(s)h^^(s) 

Gg^j  ( s )  — 

(s)g^^  (s)  (s) 
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-  transfer  function  of  the  existing  system  (jitter-to- 
position) : 


Gejit(s)- 


gi3(s) 


l+h„(s)g,i(s)f,,(s) 


-  transfer  function  of  the  existing  system  describing  the 
cross— coupling  effects  (reference— to— position  of  the  adjacent 
channel) : 

P22(s)h22(s)  gi2(s) 

G^^c(s)= - * - 

l+h22(s)g22(s)f22(s)  l+h^,  (s)  g,^  (s)  f^^  (s) 

The  following  graphs  represent  simulation  results  describing 
the  performance  of  the  described  positioning  system.  It  should 
be  noted  that  the  existing  controller  does  not  utilize  the  high 
frequency  components  of  the  position  response,  so-called  "bending 
modes",  which  presents  the  major  system  drawback. 

Figure  5-1  represents  step  and  frequency  responses  of  the 
"reference  -  position"  channel. 

Figure  5-2  shows  step  and  frequency  responses  of  the  "jitter 
-  position"  channel,  representing  system  sensitivity  to 
environmental  vibrations. 

Figure  5-3  features  step  and  frequency  responses  of  the 
"reference  -  position  of  the  coupled  channel"  thus  representing 
cross-coupling  effect. 

Figure  5-4  shows  the  control  effort  in  the  positioning 
system  while  the  system  responds  to  the  unit  step  reference 
signal. 

Figure  5-5  presents  the  family  of  step  responses  of  the 
"reference  -  position"  channel  as  the  pole  of  the  transfer 
function  of  the  controlled  plant,  normally  equal  to  -20,  varies 
from  -20  to  -10.  According  to  the  Figure,  system  performance 
deteriorates  due  to  wear  of  its  mechanical  components  thus 
exhibiting  system's  sensitivity  to  parameters  of  the  controlled 
plant. 
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3  H 


Figure  5-1 


Figxire  5-2 

STEP  EESPOHSE  OF  THE  ORIGIMAL  SVSTEM  (JITTEE  -  POSITIOW) 


uanacH-H-diQfiia  H'C'dH'c: 


Figure  5-3 
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Figure  5-5 
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6.  APPLICATION  OF  THE  MODEL  REFERENCE  APPROACH  FOR  THE 
PERFORMANCE  ENHANCEMENT  OF  A  HIGH  PRECISION  POSITION 

CONTROL  SYSTEM 


Consider  the  following  configuration  of  the  previously 
described  control  system,  modified  by  adding  a  model  reference 
(MR)  loop.  The  reference  model  symbolized  by  block  M(s) , 
represents  the  desired  behavior  of  the  entire  system.  Block  Q(s) 
is  jfesponsible  for  the  generation  of  an  additional  control 
(adaptation)  signal  W,  injected  in  the  original  control  scheme. 


Physically,  the  reference  model  is  a  computer  simulation 
code,  which  cannot  be  affected  by  any  undesirable  conditions, 
such  as  cross— coupling  between  channels,  parameter  drift,  and 
environmental  jitter.  The  output  of  the  reference  model 
therefore,  represents  the  desired  system  response  to  the  applied 
reference  signal,  reflecting  only  the  dynamics  of  the  reference 
model  and  the  reference  signal  R.  Unlike  signal  output  of 


the  control  system  X  is  dependent  upon  the  reference  signal  R, 
dynamics  of  the  "reference  -  position"  channel  of  the  system 
which,  generally  speaking,  is  different  from  the  reference  model, 
and  the  entire  set  of  "real-life"  conditions,  such  as  cross¬ 
coupling  between  channels,  parameter  drift,  and  environmental 

Therefore,  the  MR  approach  utilizes  the  feedback  control 
to  extract  the  error  signal  E=X*^^— X  and  to  convert  this 
ii^to  a  special  adaptation  signal  W=W(E)  .  When  properly 
defined,  signal  W,  injected  in  the  original  control  scheme, 
"forces"  the  output  position  X  to  follow  the  output  of  the 
reference  model.  As  with  any  system  utilizing  the  feedback 
control  principle,  the  MR  scheme  provides  a  general  purpose 
control  action  which  assures  that  error  E=x'’”-X  converges  to 
z&vo ,  or  stays  within  predictable  limits,  regardless  of  specific 
reasons  causing  the  discrepancy  between  signals  x'^”  and  X 
in  dynamics,  parameter  drift,  cross— coupling, 

jitter) . 

Introduce  transfer  matrix  M(s)  representing  the  reference 
model , 


M(s) 


m^^(s)  0  0 

0  mjj  ( s )  0 


One  can  see  that  the  reference  model  does  not  exhibit  cross- 
coupling  and  cannot  be  affected  by  environmental  jitter.  The 
particular  choice  of  transfer  functions  m^^  (s)  =m22 (s)  is  given 
below: 

m^^  (s)=1002/  (s^+lSOs+lOO^) 

Introduce  transfer  matrix  Q(s)  representing  the  adaptation 

rule. 


Q(s) 


q„(s)  0 

0  q22(s) 


While  a  general  theory  has  been  developed  for  the  synthesis  of  MR 
control  systems,  specific  adaptation  rules  can  be  successfully 
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selected  following  the  guidelines  of  [12],  and  finalized  by 
computer  simulation.  It  was  found  that  the  following  definitions 
of  transfer  functions  Q(s)  result  in  a  stable  control  system,  and 
a  good  match  between  the  outputs  of  the  reference  model  and  the 
output  of  the  control  system: 

qil(s)=q22(s)  =1000(1+.  is) /(l+.Ols) 

The  following  graphs  represent  simulation  results  describing 
the  performance  of  the  synthesized  positioning  system  with  the 
conventional  and  MR  control  loops.  It  should  be  noted  that  while 
the  conventional  controller  deals  exclusively  with  the  "main 
dynamics"  of  the  controlled  plant,  the  MR  controller  improves  the 
"main  dynamics"  of  the  system,  eliminates  "bending  modes", 
reduces  cross~coupling,  reduces  effects  of  jitter,  and  increases 
the  entire  system  robustness. 

Figure  5-6  represents  step  and  frequency  responses  of  the 
selected  reference  model.  It  should  be  noted  that  dynamic 
characteristics  of  the  reference  model  are  expected  to  be  "the 
same  or  better"  than  the  ones  of  the  original  system.  However, 
feasible  magnitude  of  the  control  effort  in  the  system  presents 
the  natural  constraint  on  the  choice  of  the  reference  model. 

Figures  5-7  shows  step  and  frequency  responses  of  the 
"reference  -  position"  channel  of  the  MR  system.  Compared  with 
the  same  characteristics  of  the  original  system  (Figure  5-8) , 
they  indicate  the  possibility  of  significant  enhancement  of  the 
system  performance,  providing  that  control  effort  (Figure  5—9)  is 
attainable  in  the  existing  system. 

Figure  5-10  shows  frequency  response  of  the  "jitter  - 
position"  channel  of  the  MR  system.  Compared  with  the  frequency 
response  of  the  original  system,  it  indicates  a  significant 
reduction  (by  20  db)  of  system  sensitivity  to  environmental 
vibrations  within  a  particular  range  of  frequencies  (10  -  10^ 
rad/ sec) .  It  is  interesting  that  unlike  adaptive  feedforward 
scheme  [2,3]  which  provides  consistent  jitter  compensation  effect 
in  spite  of  frequency  increase,  the  MR  scheme  follows  the 
feedback  principle  and,  therefore,  is  comparable  with  the 
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conventional  controller  at  high  frequencies.  At  low  frequencies, 
due  to  particular  choice  of  the  adaptation  law,  it  is  not  as 
efficient  as  the  original  feedback  controller. 

Figure  5-11  features  step  and  frequency  responses  of  the 
"reference  -  position  of  the  coupled  channel",  representing 
cross-coupling  effect  in  the  MR  system.  Compared  with  the  same 
characteristics  of  the  original  system  (Figure  5-12)  these 
results  indicate  that  the  MR  scheme,  generally  speaking,  improves 
system  performance.  However,  at  high  frequencies  frequency 
response  of  the  MR  system  is  consistently  "worse"  than  the  one  of 
the  original  system.  The  explanation  of  this  phenomenon  is  quite 
interesting.  Note  that  transfer  functions  g^2(®)  921^®) 

relate  position  outputs  of  particular  channels  to  the  control 
efforts  of  the  coupled  channels.  In  order  to  compensate  for 
"bending  modes",  the  control  effort  of  the  MR  system  contains 
high  frequency  terms,  see  Figure  5-9.  Therefore,  even  if  channel 
1  is  affected  by  a  step-type  input  signal,  transient  response  in 
channel  2  is  excited  by  the  high-frequency  control  effort  of 
channel  1.  This  results  in  the  high-frequency  resonant  peak  in 
the  frequency  response  of  the  MR  system. 

Figure  5-13  presents  the  family  of  step  responses  of  the 
"reference  -  position"  channel  in  the  MR  system  as  the  pole  of 
the  transfer  function  of  the  controlled  plant  varies  from  -20  to 
-10.  According  to  the  Figure,  application  of  the  MR  scheme 
drastically  improves  system  robustness,  making  it  virtually 
insensitive  to  parameter  drift  of  the  controlled  plant. 
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Figure  5-6 
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Figure  5-7 


STEP  RESPONSE  OF  THE  MR  SVSTEM  (REFEREMCE  -  POSITION) 


Figure  5-8 

OFIGIMAL  US  MR  SVSTEM.  STEP  RESPONSE  (REFEREMCE  -  POSITIOH) 
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Figure  5-10 
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Figure  5-12 
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Figure  5-13 
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8.  APPENDIX 


I.  The  CC  code  for  simulation,  analysis  and  design  of  the  MR 
control  system  for  500  Hz  mirror. 

echo,  Gxxl  -  main  dynamics 

k=.16  &  a=200 
Gxxl=k*a/ (s+a) 

echo,  Gxx2  -  2-nd  order  component: 
nf=3800  Sc  dr=.5  Sc  k=l 
Gxx2=k*nf^2/ (s^2+2*nf *dr*s+nf ^2) 

echo,  Gxx3  -  ocsillatory  terms 

fosc=3500  &.  dr=.03  &  k=. 000015  Sc  nf=fosc/ (l-dr^2)  ^  .  5 
gxx3  =  s*k*nf''2/  (s^2+2*nf  *dr*s+nf  ^2) 

echo,  Gxx4  -  the  drift  term 
Gxx4  = . 3 / s 

Gxx=Gxxl +Gxx2  +Gxx3  +Gxx4 
s imul , Gxx ,3, .00001,1,0, .05, 10,1, a 
freq,gxx, .1,10000,10000,2 
H=  (1+.001*S) / (1+.0001*S) 

ECHO , bode , redo 
echo, h=2+ . 02*s 
rm=l,17*6000/ (s+6000) 
gl=  (l  +  nn*h)  *  (Gxx  |  h) 
duin=Gxx*h 

g2= (l+rm*h) * (1 1 dum) 
e=nn-gxx 

simul , gl ,3, .00001,1,0, .005, 10,1, a 
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II.  The  CC  code  for  simulation,  analysis  and  design  of  the 
high  precision  MR  position  control  system  . 


echo,  g  -  transfer  function  describing  main  dynamics  (channel  1) 

g=  (s+2) / (s+20) /s 
gvl= {s+2) / (s  +  18)  /s 
gv2= (s+2) / (s+16) /s 
gv3= (s+2) / (s+14) /s 
gv4= (s+2) / {s+12) /s 
gv5= (s+2) / (s+10) /s 

echo,  b,bl  -  transfer  functions  describing  bending  modes 
echo,  in  two  adjacent  dynamic  channels 

b=100*s/ (s^2+2* .03*600*s+600^2) 
bl=126*s/ (s^2+2* . 07*800*3+800^2) 

echo,  gj  -  transfer  function  describing  jitter  effects 

gj  =310  / (s^2  +  2*.103*60*s  +  60''2)+66*s/(s^2  +  2*.  008*600*3  +  600^2) 

echo,  gc  -  transfer  function  describing  cross-coupling 
gc=2 .3*s/ ( s ^2 +2 *. 1*150*3+150^2) 

echo,  h  -  transfer  function  of  the  forward  path  filter 
h=  (s+20) / (s+40) 

echo,  f  -  transfer  function  of  the  feedback  controller 

f  =  1000/  (s+2) 

echo,  p  -  transfer  function  of  the  filter  in  the  reference 
echo,  channel 

p=1000/ (s+2) 

echo,  m  -  transfer  function  of  the  reference  model 

m=100^2/  (s"'2+2*  .  9*10  0*3  +  100^2) 

echo,  q  -  transfer  function  of  the  adaptation  mechanism 

q=1000* (l+.l*s) /(l+.01*s) 
g0= (g*h) ] f 
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echo,  transfer  function  of  the  existing  system,  ref -position : 

ge=p*gO+b 

gevl=p*g01+b 

gev2=p*g02+b 

gev3=p*g03+b 

gev4=p*g04+b 

gev5=p*g05+b 

echo,  transfer  function  of  the  existing  system,  ref-control 
echo,  effort: 
geu=p*h] (f*g) 

echo,  transfer  function  of  the  existing  system,  j itt-position : 

gejit=gj*(l| (h*g*f ) ) 

echo,  transfer  function  of  the  MR  system,  cross -coupling  effort 

gecc=geu*gc* (1 ] (h*g*f ) ) 

echo,  transfer  function  of  the  MR  system,  ref -position : 

gmr= (p*gO+m*q*gO+b) / (l+gO*q) 
gnirl=  (p*g01+m*q*g01+b)  /  (l+g01*q) 
gmr2= (p*g02+m*q*g02+b) / (l+g02*q) 
ginr3=  (p*g03 +m*q*g03 +b)  /  (l+g03*q) 
gmr4=  (p*g04+in*q*g04+b)  /  (l+g04*q) 
gmr5= (p*g05+m*q*g05+b) / (l+g05*q) 

echo,  transfer  function  of  the  MR  system,  ref-control  effort: 

g00=h| (f*g) 

gu=g00*  (p+in*q-b*q)  /  (l+gO*q) 
gul=g01* (p+m*q-bl*q) / (l+gO*q) 

echo,  transfer  function  of  the  MR  system,  j itt -position : 

gjit=gj/ (l+gO*q) 

echo,  transfer  function  of  the  MR  system,  cross-coupling  effort 

gcc=gul*gc/ (l+gO*q) 

echo,  simulation  of  a  step  response: 

simul, gecc, 3 , .001,1,0, .5, 10,1, a 

echo,  obtaining  frequency  response  and  Bode  plot: 

freq,gecc, .1,10000,5000 
bode, 3 , a 
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Figure  A-2.  Recorded  step  response  of  the  azimuth  channel 


Figure  A- 3.  Recorded  step  response  of  the  azimuth  channel 
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Figure  A-5.  Recorded  step  response  of  the  elevation  channel 
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step  response  I  step  response 


step  response  step  response 
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step  response  step  'Bsponse  step  response 


Figure  A- 11. 


Derivation  of  transfer  functions  of  the  MR  position 
control  system 


Equivalent  signal-flow  graphs  of  the  MR  system 
The  coupling  effect 


Control  effort 
of  the  coupled 
channel 


Figure  A-12.  Derivation  of  transfer  functions  or  the  MR  position 
control  system 


Equivalent  signal-flow  graphs  of  the  MR  system 
Definition  of  the  control  effort 


Derivation  of  transfer  functions  of  the  MR  position 
control  system  using  the  signal -flow  graphs 

The  ref erence- to-position  transfer  function: 

Path  #1:  P(s)  *Go(s) 

Path  #2:  M(s)  *Q  (s)  *Go  (s) 

Path  #3 :  B (s) 

Loop  #1:  -Go(s)*Q(s) 

P  (s)  *Go  (s)  +M  (s)  *Q  (s)  *Go  (s)  +B  (s) 

Gmr  (  s  )  =  - 

l+Go(s)*Q(s) 

H(s)  *G(s) 

where  Gq  ( s )  = - — - 

1+H(s)  *G(s)  *F(s) 


Let  us  check  on  the  original  signal -flow  graph: 

Path  #1:  P(s)  *H(s)  *G(s) 

Path  #2:  M (s) *Q (s) *H (s) *G (s) 

Path  #3 :  B (s) 

Loop  #1:  -H(s) *G(s) *F(s) 

Loop  #2:  -H(s)  *G(S)  *Q(s) 

P(s)  *H(s)  *G(s)  +M(s)  *Q(s)  *H(s)  *G(s)  +B(s)  [1+H(s)  *G(s)  *F(s)  ] 


1+H(s)  *G(s)  *F(s) +H(s)  *G(s)  *Q(s) 


H(s)*G(s)  H(s)*G(s) 

p(s)  -  +  M(s)Q(s)  -  +  B(s) 

l+H(s)*G(s)*F(s)  l+H(s)*G(s)*F(s) 

1+H  (s)  *G  (s)  *F  (s)  H(s)*G(s) 

-  +  Q(s)  - 

1+H  (s)  *G(s)*F(s)  l+H(s)*G(s)*F(s) 


Definition  of  the  control  effort: 


Path  #1:  P(s) *H(s) 

Path  #2:  M(s)  *Q(s)  *H(s) 

Path  #3 :  -B (s) *Q(s) *H{s) 

Loop  #1:  -H(s) *G(s) *F(s) 

Loop  #2  :  -H(s)  *G(s)  *Q(s) 

P(s)  *H(s)  +M(s)  *Q(s)  *H(s)  -B(s)  *Q(s)  *H(s) 

Gjj  { s )  = 

1+H(s)  *G(s)  *F(s)  +H(s)  *G(s)  *Q{s) 
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H(s) 

introduce  Gi(s)  = - - -  ,  then 

1+H(s)  *G(s)  *F(s) 

Note  that  G]_  (s)  G  (s)  =Gq  (s)  ,  therefore, 


Gu{s)=G,(s) 


P  (s)  +M  (s)  *Q  (s)  -B  (s)  *Q  (s) 
1+Q  (s) *Go  (s) 


Definition  of  the  jitter  effects: 

Path  #1 :  Gj  (s) 

Loop  #1 :  -Go  (s) *Q (s) 

Gj  (s) 

Gjit  (  S  )  =  - 

l+Go(s)  *Q(s) 


Analysis  of  the  cross-coupling  effect 

Path  #1:  G^(s) 

Loop  #1  ;  -Go  (s) *Q (s) 

G,  (s) 

Gqc  ( s )  =Gu  ( s )  - - — 

l+Go(s)  *Q(s) 


»U.S.  GOVERNMENT  PRJNTING  OFFrCE:  1997-509-127-61051 
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MISSION 

OF 

ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 

a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Material 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence, 
reliability  science,  electro-magnetic  technology,  photonics,  signal 
processing,  and  computational  science. 

The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


